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ABSTRACT 


Corrugated  metal  (and  "doubly  corrugated  metal")  la 
an  anisotropic  boundary  which  supports  only  a hybrid  surface 
wave  mode  in  directions  other  than  parallel  or  perpendicular 
to  the  corrugations.  The  mode  phase  velocity  is  a function 
of  direction.  To  determine  the  excitation  efficiency  of  these 
modes  a three  dimensional  analysis  is  performed  generalizing 
the  results  of  Cullen  [l 1 and  utilizing  the  author's  previous 

work'2]-  v kj 

The  surface  waves  are  shown  to  possess  many  of  the 
properties  of  plane  waves  in  a two  dimepsional  anisotropic 
medium  - notably  that  the  energy  propagation  is  radial  and 
not  generally  normal  to  the  wave  fropxs.  Radial  propagation, 
however,  does  not  Imply  that  the  primary  pattern  of  a feed  and 
the  resulting  far  pattern  of  the  surface  wave  are  the  same. 

The  hybrid  plane  wave  component  off  the  spectrum  of  the  source 
which  propagates  in  a direction  yv-parallel  to  the  ground  plane 
with  the  natural  surface  wave  phase  velocity  of  the  boundary 
in  this  direction  is  observed  in  the  far  field  in  a direction 
[0  not  generally  equal  to  Y i Moreover  the  excitation  efficiency 
is  not  independent  of  VC.1l  Therefore,  the  possibility  of  focussing 
by  an  anisotrogifi_JioundaryJ  remains . 

ulk  f< 


Part  II  Cullen's  formula  for  surface  wave  excitation 


efficiency  over  a dielectric  slab  and  ground  plane  is  also 
generalized  to  thick  slabs.  The  wide  range  of  efficiencies 
obtainable  as  a function  of  slab  thickness,  dielectric  constant, 
and  source  height  is  related  to  many  of  the  difficulties  ex* 
perienced  in  constructing  surface  wave  antennas. 
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PART  I 

EXCITATION  OF  SURFACE  WAVES  IN  FREE  SPACE 
OVER  AN  INFINITE  PLANE  ANISOTROPIC  BOUNDARY 


rionucTioN 


Cullen1-  J has  considered  Che  excitation  efficiency  of  a 
surface  wave  over  an  infinite  plane  corrugated  metal  boundary 
launched  by  an  infinite  line  of  magnetic  current  parallel  to  the 
corrugations.  In  this  paper  we  consider  the  more  general  problem: 
of  an  arbitrary  source  over  corrugated  metal  and  other  boundaries 
to  be  described. 

From  a three  dimensional  point  of  view  corrugated  metal  is  an 
anisotropic  boundary.  The  surface  waves  excited  on  it  have  phase 
velocity  dependent  on  the  direction  of  propagation  and  many  of  the 
other  properties  of  plane  waves  in  a two  dimensional  anisotropic 
medium.  The  surface  wave  front  created  by  a point  source  near  an 
anisotropic  boundary  is,  in  fact,  a more  complicated  curve  than  the 
elliptical  wave  front  of  a point  source  in  a two  dimensional  aniso- 
tropic medium.  The  boundary  anisotropy  may  cause  focussing  or  de- 
focussing of  the  source  - a phenomenon  with  no  counterpart  in  Cullen's 
two  dimensional  problem. 

In  practical  microwave  work  one  is  sometimes  confronted  with 
a three  dimensional  problem  in  which  there  is  a plane  of  symmetry 
both  with  respect  to  the  geometry  and  with  respect  to  the  incident 
field.  In  this  case  it  is  common  to  assume  that  the  solution 
to  the  problem  in  the  plane  of  symmetry  is  not  very  different 
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from  that  of  the  "two'  dimens ionalized"  problem  which  one  obtains 
by  assuming  that  the  incident  field  and  the  geometry  are  the 
same  in  the  plane  of  symmetry  as  in  the  original  problem  and 
are  independent  of  the  coordinate  perpendicular  to  this  plane. 
This  assumption  is  implicit  in  [ 1 ] , for  although  in  the  main 
body  of  the  paper  Cullen  clearly  states  that  he  is  solving  a 
two  dimensional  problem,  in  his  section  on  experimental  veri- 
fication he  describes  a three  dimensional  set  up.  His  ex- 
perimental source  is  a half  wave  long  slot  rather  than  the 
infinite  magnetic  line  current  assumed  in  his  theory.  The 
failure  to  distinguish  between  these  two  sources,  might  very 
well  not  offend  the  microwave  intuition  of  anyone.  How  - 
ever,  in  this  paper  the  efficiencies  of  surface  wave  excita- 
tion are  compared  for  the  two  sources  and  they  are  found  to 
differ  both  because  of  Inherent  differences  between  two  and 
three  dimensions  (chis  applies  both  to  the  corrugated  and  to 
the  dielectric  clad  ground  plane)  and  also  in  the  case  of  the 
corrugated  ground  plane,  because  of  the  anisotropy.  In  certain 
cases  the  difference  may  be  numerically  appreciable.* 

If  the  corrugations  are  parallel  to  the  plane  x » 0 
and  the  z axis  extends  out  of  the  grooves  then  in  the  two 
dimensional  case  the  boundary  is  characterized  by  a normal 
impedance  condition: 


(1) 


z-o 


-Z 


I 


* See  Figures  9 a,  b,  c,  and  d. 
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- **  i 

If  the  thickness  of  the  teeth  of  the  corrugations  is  negligible 
compared  to  the  gap  between  them  then  ] 


(2) 


Z 


tan(kd) 


where  d is  the  depth  of  the  grooves  and  k « a>  is  the 

free  space  propagation  constant^.  The  same  elementary  con- 
siderations which  lead  to  (2)  show  that 


Before  proceeding  it  is  useful  to  generalize  the  problem 
from  the  case  of  ordinary  or  single  corrugations  to  the  case  of  a 
"doubly  corrugated  surface"  shown  in  Figure  1. 
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Such  a surface  may  be  made  by  milling  on  a flat  metal  plate 
a set  of  parallel  grooves  of  depth  d , width  g , and  tooth 
width  t , then  rotating  the  plate  90°  and  milling  a second 
set  of  grooves  of  depth  0,  width  G,  and  tooth  width  T.  If 
the  x direction  is  taken  parallel  to  the  first  set  and  the  y 
direction  parallel  to  the  second  set,  then  to  replace  (1)  - (3) 
we  have  the  following  boundary  conditions 

(4)  E . * 

ir  2 , o - -z  - ^ pF 

(5)  h*  - tan(kO)  . 

Here  G,  T,  g,  and  t must  all  be  assumed  small  compared  to 
the  wavelength.  Obviously,  if  we  take  D«  0 or  G * 0,  we 
return  to  the  case  of  single  corrugations. 

SOLUTION  OF  A PLANE  WAVE  INCIDENT  ON  A DOUBLY  CORRUGATED  BOUNDARY 
The  two  normal  impedance  conditions  (4)  and  (5) 
determine  the  scattered  field  when  any  field  is  incident 
upon  the  surface.  This  will  be  shown  in  general  later.  Let 
us  first  assume  that  a TM  plane  wave  is  incident  at  angle  0 
to  the  outward  surface  normal  with  the  plane  of  incidence  at 
angle  V to  the  positive  x axis.  Let 

(6)  p « ksinO,  w ■ kcosd  **  J k^-p^ 

and  let  ^ , 0,  and  z be  the  cylindrical  coordinate  system 
associated  with  the  rectangular  coordinates  in  the  usual  way. 

Then  the  incident  plane  wave  field  may  be  written  as 
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(7) 


H 


e-iwz+ijy»cos  (0-  V*) 
w” 


; ' <i 

(sinV,  -cos V,  0)  '] 


«t  - e-i'“+iWcos«-V'><cosV,  sin  V,  -g-) 

We  use  here  the  notation  (a,  b,  c)  for  the  x,  y,  and  z 

components  respectively  of  a vector. 

In  order  to  satisfy  the  boundary  conditions  (4)  and 

(5)  it  is  necessary  in  general  to  postulate  two  plane  waves 

reflected  from  the  surface,  one  TE  and  one  TM,  rather  than 

simply  a TM  reflected  plane  wave  as  is  sufficient  in  the  case 

of  an  isotropic  boundary  condition.  Mathematically,  this  is 

true  simply  because  we  need  in  general  two  reflection  coefficients 

and  R2  to  satisfy  the  two  conditions  (4)  and  (5).  Let  the 

reflected  TM  wave  be 

R giwz+ip^cos  (0-  Y ) 

(8)  H - -i-  (sin  V,  -cosV,  0) 

w 


1 


i 


<ocE  - R1etwz+1PA°‘(«-V)(.cosV>  _slny  _E_) 
and  the  reflected  TE  wave  be 

(9)  kH  - R2eiw*+1Wco,<*-V>(co.y,  sin  Y,  - -£-) 


Js-e  . h ,iw+iiy«>.«-Y)(,in y,  ,cosV(  0). 

Each  of  these  waves  satisfies  Maxwell's  aquations  and  propagates 
outward  to  infinity.  The  total  field  will  satisfy  the  boundary 
conditions  if,  as  we  find  upon  substitution  in  (4)~(5), 
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(10)  -Rj(-^j-  + + “§(1  + 5^~)cos  V 


sin  V <-±-  - s^-) 

o 


MHr + sb)c0S  v + “c2-  (1  + Si)5ln  v 


C0,v  <sk-^3>  • 


This  pair  of  simultaneous  linear  equations  may 
readily  be  solved  to  yield 

(11)  (r  +X0)(*X+l)sin2Y  + ( Y + X)(yXQ+l)cosV 

, « 

(y-xo)(yx+i)sin2v  + ( y - x)(/xo+i)cos2v 

2(XQsin2v  + X cos2  y + V XXQ) 
y^(X  cos2y+  Xsin2 V ) + y(l-XX  )-Xcos^y-X  sin2y 


21  sin  Ycos  V y (X  -X) 

C12)  IL>  * o 

^ < Y -XQ) ( y X+l)sin2 v + ( y -X) ( y xo+i)cos2  Y 


where  we  have  set 


(13) 


y-iS- 


* « i 


/?*• 


X - iZ 
o o 


In  the  lossless,  doubly  corrugated  surface  case,  X and  XQ 
are  real.  If  V is  real  and  positive  the  reflected  waves 
are  surface  waves  traveling  in  direction Y • 

We  have  thus  solved  the  problem  of  reflection  of 
an  Incident  TM  plane  wave.  In  order  to  do  the  same  for  a 


TE  plane  wave,  we  can  use  the  duality  principal  on  all  of 
the  preceeding  equations.  The  duality  principal  states  that 
these  equations  remain  correct  under  the  substitutions: 

(15)  E »H 


if  these  substitutions  are  made  wherever  these  symbols  occur. 
The  results  give  the  solution  for  the  TE  plane  wave.  An 
arbitrary  source  can  be  resolved  into  two  sources  producing 
a continuum  of  TE  and  1M  plane  waves  respectively.  This  is 
proved  in  [2].  We  will  next  show  how  to  find  the  total 
solution  for  the  TM  part  of  a general  source.  When  this  has 
been  done,  we  can  apply  the  substitutions  (15)  to  the  TM 
solution  in  order  to  yield  the  solution  for  the  TE  part  of 
the  source,  and  thus  the  total  solution.  Bearing  in  mind  that 
there  is  no  essential  loss  in  generality  thereby,  we  confine 
our  attention  exclusively  to  a TM  incident  field  for  the 
remainder  of  the  paper.  Before  doing  so,  however,  it  is 
worth  noting  that  the  poles  of  and  R2  as  functions  of  p 
are  unchanged  by  the  substitutions  (15).  This  implies,  as 
we  shall  see,  that  in  general  the  surface  waves  which  can  be 
produced  by  TE  and  TM  sources  are  the  same.  In  any  case  no 
surface  waves  exist  other  than  those  corresponding  to  the 
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poles  of  R^  (or  R^  which  has  the  same  poles  as  R^). 

SURFACE  HAVES 

When  the  denominator  of  in  (11)  is  zero  a free 
solution  to  Maxwell's  equations  can  exist  which  satisfies 
the  boundary  conditions  and  has  many  of  the  features  usually 
associated  with  a surface  wave.  We  shall  derive  these 
surface  wave  fields  in  a later  section.  Let  us,  for  the 
present,  however,  simply  examine  when  they  can  occur.  The 
condition  for  a surface  wave  is  that  its  vertical  attenuation 
constant  shall  be  YQ,  the  particular  value  of  Y which  is 
positive  and  satisfies: 

(16a)  > (X  cos2y+XsinV)+  V0(l-XX)-XcosV -X  sin2y-  - 0 


1+iTX, 


(16b)  tan2y  « f(  V )g(  Yq),  £(Y)-2-^JL  ,g(/) 

° ° 1 + if  X • 

(16a)  and  (16b)  are  quadratic  equations  for  Y Q,  having  at 
most  two  roots.  These  are: 


(17) 


XX  -1  - /(1+XXJ2  + (X-X  )2sin22Yf 


2 I X + (XQ  - X)cos4 y 


When  the  denominator  vanishes,  tan = -yP—  and 
„ X + X (l-tan2y) 

(18)  Y m ——————  a — - ■■  - — m 

° (X-XQ)(1-XX0)  (1+Xztanzy)  (1+tan2  y) 

If  X and  XQ  are  both  positive,  then  the  only  positive  value 
of  Y corresponds  to  the  positive  sign  for  the  radical.  If 


X and  XQ  are  both  negative,  then  the  negative  sign  is  required 
for  the  radical.  If > one  is  positive  and  the  other  negative, 
either  sign  may  be  required.  In  order  to  determine  the 
number  of  positive  solutions  of  (17)  and  some  of  their 
properties  as  functions  of  V we  now  consider  various  cases 
separately. 

(i)  Pure  Inductive  Case 

Assume  X and  XQ  are  both  positive.  Without  loss  in 
generality,  let  X^XQ.  For  each  V , there  is  then  one  and 
only  one  Y satisfying  (16).  We  may  see  this  as  follows: 

The  left  member  of  (16b)  is  always  positive.  The  right 
member  can  only  be  positive  if  YQ  is  between  XQ  and  X.  In 
the  range  XQ Y ^ X,  both  f(Y)  and  g (/  ) are  positive  de- 
creasing functions.  Therefore,  their  product  is  also  a mono- 
tone function  which  decreases  from  infinity  to  0 in  this  interval 
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Since  the  left  member  is  a monotone  increasing  function 


for  O^y^r/2,  there  can  be  only  one  solution 


Furthermore,  this  solution  equals  X when  V ■ 0 and  decreases 
continuously,  as  Y increases,  reaching  a value  XQ  when 
V"  “ *72.  The  values  repeat  periodically  in  the  other 
quadrants  - as  indeed  they  must  - for,  from  the  symmetry 
of  the  problem,  there  is  no  distinction  between  any  angle 
Y and  the  related  angles  - Y , t r-  Y , or  i r + Y . There- 


fore yQ  is  a maximum  when  V * 0 and  ir  and  a minimum  when 


Y * v/2  and  3ir/2,  and  accordingly 

when  0 Z w/2 


, N . ay0<Y> 

(19)  « S-jA— 


0 when  V ■ 0,  ir/2,  r,  or  3r/2. 


(ii)  Mixed  Case;  AikO,  XQ  j£L  0 


(iia)  XXQ  ^ -1  or  X*. — 


In  this  case,  in  the  interval  0 both  f ( y ) 

and  g(y)  are  positive  decreasing  functions,  hence  so  is 
their  product.  Equation  (16)  has  a unique  root  here  if  and 
only  if 


(20) 


tan  V"  ~ 


For  this  root  y*  0.  In  the  interval o-  £ t & X, 

° Ao 

f ( / ) is  positive  and  g(Y)  is  negative,  so  that  there  are 

no  solutions  of  (16).  If  )f  ^ X,  both  f(Y)  and  g(  y)  are 

negative  decreasing  functions  so  that  their  product  is  a 

positive  increasing  function.  Equation  (16)  then  has  a 

unique  root  in  this  interval  if  and  only  if 

^ 7 x« 

\ A ^ 4- ^ ^ O 


(21) 


0 ^ tanz  - 


and  for  this  root  y’  ^ 

o 

as  follows:  Let  m « min 


0.  These  facts  may  be  summarized 


f-v  -xl 


Then  m 


We  conclude  that  if  tan  V"  £ m,  there  are  two  values  of 
satisfying  (16),  one  in  the  interval  0 y -1/XQ 

and  one  satisfying  X 4s  Y . If  m ^ tan2  V ^ there 

is  only  one  value  of  Y_  satisfying  (16)  and  it  lies  in 
the  interval  0 Y - 1/XQ  if  m ■ - -y  , or  in 


y £ X if  m - - 


- 1/XQ  if  m j—  , or  in 

If  tan2Y^,-~  , (16)  has  no 


solutions.  The  root  in  the  interval  0 4^  - 1/XQ  is 

a decreasing  function  of  Y :,  i.e.,  y' ^ 0,  and  the 

o 

root  in  the  interval  X £ y is  an  increasing  function  of 

Y : i.e.,  ^ 0 . This  case  is  illustrated  in  Figure  7 

o 

with  X - 1,  X0  - -4. 
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(iib)  XXG  ^ -1 


The  conclusions  in  this  case  are  similar  to  (iia). 

If  tan2y  £ m,  there  are  two  solutions  to  (16),  one  in 
the  interval  0 ^ Y £ X and  one  in  the  interval  -1/XQ  £ / 
If  o ^ tan2  V £ - , there  is  only  one  solution.  It 


satisfies  -1/XQ  * y if  m - - 


or  0 <£  Y <£  X if 


m . Both  cases  (ii  a)  and  (ii  b)  may  be  summarised 

together  if  we  let  M ' - max  ^X,  - > m “ min 

and  m'  - min  £x,  -1/XQ  } . Then  if  tan2  y £ m,  there  are 
two  solutions  to  (16),  one  in  the  interval  0£  Y £ m* 

(for  which  y'Q  0)  and  one  in  the  interval  M'  *£  V (for 

which  ro  ^ 0)-  If  « ^ tan2  -J-  * there  is  only 
one  solution.  It  satisfies  M'^s  Y if  m ■ — v or 


0 ^ y^a*  If  b ■ - 


If  tan2y  , there 


ere  no  solutions  to  (16). 

i 

(iii)  Mixed  Case:  XQ  ^ 0,  X zl  0. 

Here  if  we  let  M"  ■ max  ^XQ, , 

“ min  £ *0> X“J  • Then  if  * tan2  y , there  are 

two  distinct  solutions  to  (16),  one  in  the  interval 

0 y £=.  m"  (for  which  ^ 0)  and  one  in  the  interval 

M Y (for  which  x'Q  ^ 0) . If  m ^ tan2  y-  * , there 

is  one  solution.  It  satisfies  M"  4 y>if  m , or 

X ° 

0 ^ Y »”>  if  m ■ - — jj—  . If  tan2  y m there  are  no 
solutions . 

(iv)  Pure  Capacitive  Case:  X 0,  X 0 
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In  this  case,  for  each  V there  is  a unique  positive  YQ 
satisfying  (16).  This  value  of  YQ  lies  between  -1/X  and 
-1/XQ.  If  we  assume  XQ  ^ X,  then  jL  0 . 

FIELD  PRODUCED  BY  A SOURCE  OVER  A DOUBLY  CORRUGATED  BOUNDARY 


Equation  (38)  of  reference  [2]  gives  the  represen- 
tation of  a more  or  less  arbitrary  incident  TE  field  as  a 
superposition  of  plane  and  hybrid  waves  weighted  according 
to  the  far  field  of  the  0-component  of  E for  the  source  in 
free  space.  Equation  (125)  of  [2]  gives  a similar  repre- 
sentation for  the  scattered  field.  In  order  to  utilize 
these  results  here  we  must  first  employ  the  duality  sub- 
stitutions (15).  We  must  then  make  a change  in  coordinates 
(in  [2]  positive  z points  into  the  boundary).  Finally  we 
must  Include  the  pair  of  reflected  plane  waves  instead  of 
the  single  plane  wave  reflected  from  an  isotropic  boundary. 
If  these  steps  are  performed  correctly,  we  obtain  the  total 
field,  incident  plus  scattered,  for  a TM  source: 

(22)  " 2t  ip/»cos(0-V)+iwz 


(— ■ 

I H 


dpd  V^ks  inV’(U++R1U J 


-R^wcosyU_,  kcosV(U++RjU__)  - R2WS inY U_,  R^pU^  ^ 
• ow  tpAioa  (0-V)+iwz 

4 / *»  n • t 


E.  i ( Ee 

2tu>€  j ) w 


dpdvjwcos  V (U++RXUJ 


-kRjSinV U_,  wsin V (U++RXU_)  + kR2cosV  U_  , 


-p(U+  + RXUJ  j 


TCruWifAl  PCCCAD 
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The  symbolism  U+  and  U_  has  the  following  meaning:  the 

0 component  of  H , when  the  source  is  in  free  space, 

« 

assumes  the  form 

_ikr  5 

(23)  F(9,0>S—  + 0(-4-)  as  r — >oo 

r r* 

where  F(9,0),  the  "far  field",  may  be  rewritten  in  terms 

of  p « ksinO  and  V « 0 as  a new  function 

u(PiV)  “ F(sin“^p/k,  V )•  An  ambiguity  exists  in  this 

definition  since  the  inverse  sine  has  two  values: 

ir*0  as  well  as  9.  The  far  fields  in  these  two  directions 


are  not  generally  the  same,  of  course.  (An  exception  in 
the  plane  0 « 0 is  the  magnetic  line  current  along  the 
y axis.  Here  F(0,O)  is  a constant).  U(p,  y ) is  thus 
generally  a double  valued  function.  In  (22)  we  have 
let  U+  stand  for  the  value  of  U(p, y ) corresponding  to 
directions  for  which  9 is  less  than  v/2  and  U_,  the 
other  value,  for  which  9 exceeds  r/2.  Our  problem  now 
is  to  obtain  both  the  radiation  field  and  the  surface 
wave  field  from  (22).  The  former  is  the  asymptotic 
value  of  these  integrals  for  large  r,  where  />  * rsintf, 
z - rcosO.  The  latter  is  the  asymptotic  value  when  is 
large  and  z is  fixed.  Faced  with  the  same  problem  for 
the  iso tropic  boundary  case,  the  author  required  con- 
siderable mathematical  analysis  in  reference  [2J  to 
obtain  just  the  surface  wave  field  rigorously.  The 
method  employed  there  cannot  apply  directly  here  since 
now  the  reflection  coefficients  have  an  angular 


16. 


dependence:  and  are  functions  of  V* 

Instead  here  we  outline  a proof  which  it  is  believed 

« 

could  be  made  rigorous  though  not  without  great  effort.  We 
present  it  simply  as  a plausible  argument.  We  begin  by  writing 
each  of  the  Integrals  of  (22)  as  a sum  of  two  terms , one  contain- 
ing the  V Integration  from  -t/2  + 0 to  v/2  + 0,  and  the  other  the 
integration  over  the  remaining  interval  from  r/2  + 0 to  3v/2  + 0. 

We  Indicate  this  symbolically  with  H as 

? 0+t/2  • 0+3»/2 

(24)  H - Hr  + H2,  Hj  - ) { , h2  - ) j 

0-r/2  * i+w/2 


since  there  is  no  need  here  to  bore  the  typist  with  the  lengthy 
integrands.  Now  we  observe  that  U(p,Y)  - U(-p,  Y + *)»  since 
the  direction  (0,0)  is  the  same  as  (-0,0+r).  Thus,  if  in  the 
Integrand  of  H2  we  make  the  change  of  variable  v ' ■ V + v. 


p*  a -p. 

we  find  that 

! 0+t/2 

(25) 

H2-J  J 

0-f/2 

to  that 

r 0 + w 2 

(26) 

■- ) ] 

— 0-TT/2 

(The  details  of  Indentation  of  the  p integration 
contour  about  the  poles  of  R^  and  R2,  of  course,  must  be  properly 
hendled).  In  the  p integration  we  have  succeeded,  therefore. 


T 
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in  obtaining  an  Infinite  contour  which  may  - if  the  Integrand 
behave#  properly  at  infinity  - be  deformed  around  the  poles  and 
branch  cut.  In  the  V integration,  we  have  avoided  the  stationary 
phase  point  in  the  backward  direction  V - 0 + r and  are  left 
with  the  only  physically  reasonable  stationary  phase  point  at 
Y An  estimate  of  U(p,  y ) for  large  p is  required.  It 

is  not  hard  to  show,  from  the  formulas  of  [2]  , that 


(27) 


U(p,V) 


/o  Jm(p)  - f^e(p)/r 
0(e  0 - 1 - ),  p — *00 


where  /°Q  is  the  radius  of  the  smallest  sphere  containing  the 
source  and  S is  the  minimum  distance  of  the  source  to  the  plane 
z ■ 0.  A formula  similar  to  (27)  would  be  required  in  a rigorous 
proof. 

RADIATION  FIELD 

The  radiation  field  follows  from  a stationary  phase 
evaluation  of  the  integrals.  When  the  observation  direction 
is  (0,0)  the  stationary  phase  point  occurs  when  V * 0# 

P - ksinO  (and  w - +kcos9).  The  results  may  be  most  simply 
expressed  in  terms  of  the  0 and  0 components  of  the  field  as 
follows : 

(28)  H#  - F(9,0)  + R]F(r-9,0) 

- -*^(*-0,0) 

*f  - -f*-  KjFtr-M) 

*•  * * *!?<*-•,#» 


aaAlifl 


i 
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where  end  Rj  are  evaluated  at  / - -icosO,  V « 0: 

(29)  Rj  - ^coa29  ~ XQ)sin2  0 + (Xoco82O-X)cos20  + icos0(l+XXo) 
(icos0+Xo)(l-iXcos0)sin20  + (icos0+X)(l-iXocosO)cos20 

^ 2(iXXQcos0  - XQsin20  - Xcos20) 

(Xocos20  + Xsin20)cos20  - icos0(XXo-l)  + Xco820+XQ8in20 

2sln0cos0cos9 (X-XQ) 


(Xqcos20  + Xsin20)co829  - icos0(XXo-l)  + Xcos20  + XQsin20 
The  power  flow  between  two  meridians:  0 Z.  0Z.  0,  is 

Cr/2  *1  °.  N1 

(30)  Pr  - — j ( sinOdO  f d0(ExH)  • n i 

0 ' 

Fo 

where  n is  the  unit  radial  vector.  Substitution  of  (28) 
yields 

(31)  W^2  ^l 

Fr  " “J -J^F~  \ sln9d0  \ d?J  | F(0.*)+R1F(»'-8,0)|2  + 

0 *.  I 

jRjFCr-e.l/j  . 

The  first  term  in  the  integrand  corresponds  to  the  unde  — 
polarized  field,  the  second  term  to  the  depolarized  field. 

The  two  fields  are  orthogonal  and  their  energies  are  additive. 
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We  observe  that  the  poles  of  and  R2  are  functions 
of  V • Hence  we  perform  the  p integration  first.  While, 
as  mentioned  above,  our  derivation  is  not  meant  to  be  completely 
rigorous,  it  is  worthwhile  to  determine  carefully  whether  the 
poles  are  actually  crossed  when  the  p integration  contour  is 
deformed.  The  whole  question  of  the  existence  of  the  surface 
waves  is  at  stake  here.  For  this  purpose,  as  is  customary  in 
similar  problems,  let  us  introduce  a little  conductivity  into 
the  medium.  That  is,  we  assume  that*£n(k)  (the  imaginary  part 
of  k)  is  slightly  positive.  The  integrands  of  (22)  have  branch 
points  at  p ■ k and  -k.  The  integration  contour  is  along  the 
real  p axis  and  we  have  taken  branch  cuts  as  shown  in  the 


figure  (the  same  as  in  ref.  [ 2]  ). 

vMp)  l 


pole 


i^e(p) 


branch  cuts 


FICURE  NO.  2 


He  must  deform  the  contour  into  the  upper  half  plane  in  order 
to  safely  neglect  the  contribution  from  infinity.  He  shall  obtain 
a residue  from  a pole,  therefore,  only  if  it  lies  in  the  upper  half 
plane.  Let  us  show  that  all  the  poles  which  are  real  when«lm(k)  * 0 
do  lie  in  the  upper  half  plane  when  tlm(k)  is  positive  and  approach 
the  real  axis  as«Ai(k)  approaches  zero.  Equations  (4),  (5),  and  (13) 
imply  that  both  X and  XQ  take  on  a slight  positive  imaginary  part  when 
k does.  Suppose  for  the  moment  that  V is  fixed,  then  (16)  defines 
YQ  as  a function  of  k . When  A(k)  is  zero,  of  course,  YQ  is 
real.  Implicit  differentiation  of  (16)  with  respect  to «/m(k)  (bear- 
ing in  mind  that  all  three  quantities  VQ>  X,  and  XQ  are  functions 
of  Im(k))  shows  that  the  derivation  of  J^m(  yQ)  with  respect  to 
c/n(k)  is  always  positive  when  YQ,  X,  and  XQ  are  real.  Thus  the 
introduction  of  slight  conductivity  moves  the  pole  (or  poles)  in  the 
/ -plane  from  the  real  axis  slightly  into  the  upper  half  plane. 

The  cut  p plane  of  Figure  2 is  mapped  into  the  region  lying  to 
the  right  of  the  Image  of  the  branch  cuts  in  the  Y plane  (Fig.  3), 


■ -H* 
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Since  p - k >/l  + )r  , 


(32) 


^ tPi n(k) 


'1  + * 


} tPm i(k) 


By  setting  elm(k)  equal  to  zero  in  (32),  if  YQ  is  positive,  we 
see  that  the  left  member  is  positive.  Every  pole  in  the 
plane  which  maps  into  the  cut  p plane  corresponds  to  positive 
and  hence  maps  into  the  upper  half  of  the  p plane.  Thus 
the  residues  of  all  poles  are  actually  obtained. 

Taking  the  limit  as  tAn(k)  vanishes,  we  find  from  the 
residue  theorem  the  surface  wave  terms 


0+ir/2  , - 

(33)  H - l f ^2  B(k/f  ,V)e‘k  ^oz+ikAtfocos (-r  slnV 
»-w/2  r°  .1 

-i  v2co8Y > r^coeY  - i Y0rzsinY,  r2/tfQ)dY 

0+w/2 

Es  - £ ( U(k/^o.V)e-kyoz+iV,^oco8(0"^  >(i  YorlCos  V 

0-*/2 

-r2sinV,  i Y^sinV  + r2cosV,  ■r1/0)dV> 

where,  we  have  used  our  previously  defined  vector  notation  and 
have  let  ^ be  the  "equivalent  index  of  refraction"  of  the 
surface  wave: 


(34) 


A -A  + i. 


(35) 


resRi 


P-k/o 


2k  /0(Xcos‘iv'+X0sinzy+  yQXX0) 

(2  yo(Xsin2y+X0co»2v)  - XXQ+l)  /l+  /D2 


resR, 


2isiny’cosV’k  y *(X  -X) 


(2  yo(Xsin2y+  Xocos2y  ) - XXQ  + £ /i  7T2 


If  there  is  mere  than  one  positive  value  of  Yq  satisfying  (16) 

(as  there  may  be  in  the  mixed  case  XXQ^.0),  then  there  is  one 
term  of  the  form  (33)  for  each  value  of  Y . 

We  next  evaluate  the  integrals  for  large /°  by  the  method 
of  stationary  phase.  It  is  at  this  point  that  an  important 
difference  enters  between  the  present  problem  on  the  one  hand, 
and  the  two  dimensional  or  Isotropic  boundary  problems  on  the 
other  hand.  The  difference  is  that  in  the  exponential  coefficient 
k^coi^'V)  of  , the  factor</^Q  depends  on  Y • The  stationary 
phase  point,  thus,  in  general  does  not  occur  at  V * 0.  It 


occurs  at  the  root  of 


+ 1 cos 


If  we  call  this  root  V (0)  Chen,  by  carrying  out  the 


differentiation  we  find 


0 - Y - tan' 


•i/liais 
u + K 


Y - tan' 


ft). 


. ‘s&u&m 


where  V0  * ) #0/<)  V*  • and  c^e  principal  value  of  the  inverse 
tangent  is  implied  in  (37) . Since  Y is  given  as  a function 
of  Y by  (17) , for  each  0 the  stationary  phase  point  y may  be 
found,  in  principal,  by  means  of  (37).  Calculations  for  simple 
corrugations  (XQ  « 0)  are  graphed  in  Figure  4.  It  appears 
possible  in  certain  cases,  to  have  more  than  one  stationary 
phase  point,  i.e.,  two  or  more  values  ofV  for  a given  0.  We 
shall  consider  first  those  cases  where  only  a single  stationary 
phase  point  occurs. 

For  the  stationary  phase  evaluation  we  also  require 


(38) 


s<0) 


yo 


I 


Y-YW 


'2  2 2 ,2  2 '2 

K - y,)-2  K I 


[a+  v*><(i+  >l>2  + Vo  *o2)] 1/2 


We  obtain  finally  the  surface  wave  field  as 


(1+i)  i-T—  . y kz+il^  i 008(0-^)^ 

(39)  H.  * y0  JWi  U(k/Vy)e  (-r^inV 


-i^0r2cosy,  r^osy-i  ^rjSinV,  r2/<f>0) 


U*U/*o* 
0»€  * 


U(k/dD,  y ) e"  <>kz+ik/^oco8  * (i  J^r^cosy 


-r2siny,  i y^siny  + r2cosy,  - r1/d,Q) 


* This  is  shown  in  Figure  7 for  X ■ 


Sity) 


cosi<p-V) 


L 
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FIGURE  4a 

INDUCTIVE  CASE:  X - 1,  X{  - 0.  Showing  Far 
Field  Wave  Front  and  Relation  of  Initial  Energy 
Flow  Direction  V and  Far  Field  Energy  Flow  0. 
(Isotropic  Phase  Source) 

, TECHNICAL,  BESCABCW  GROUP 


^ ♦ fo  cos  (y-yr) 

FIGURE  4c 

INDUCTIVE  CASE:  X - 5,  XQ  - 0.  Showing  Far 
Field  Wave  FronC  and  Relation  of  Initial  Energy 
Flow  Direction  V and  Far  Field  Energy  Flow  0. 
(Xeotroplc  Phase  Source) 
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The  power  flow  between  the  meridians  0QZL  0^0^  is 
/p  r OO  0X  \ 

(40)  dz  j />d0(Eo  X H„)  $ r 


/*d«(Es  x 


5s>j  1 


lx  3 / 

3Kc  f jfffs  | U(k/ , V)/  2d«  kri2-r|)co6  Y -21 

(r^-rjJslnV  + 21  J^r^cos  y,  o) 

where  V is  defined  as  a function  of  0 by  (37).  The  surface 
wave  power  flow  is  thus  parallel  to  the  surface.  The  radial 
component  over  an  element  d0  is 


<4i)p  _ w/fl lu(yy,y)|: 

sf  2ko>e  jsj  yJ 


<r  1”  r2*cos  V )+2i*orjr2sin  (0-y  ) 


The  azimuthal  or  0-component  is 


<42V  I^///Uyy)/  ^ (r22-r12)»ln«-y)t21>ro»1r2co.(0-V^ 

By  direct  but  lengthy  computation  (details  in  Appendix  I),  we 
can  prove  that 

(43)  (r22-r12)sin(0-y)  + 2i>’or1r2cos(0-y )»  0, 


As  s consequence. 
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1 

f 

t 


j 

i 


(44) 

and 

(45) 


P.,-0 

, _p  M„3lu(M,y)|2<ri2-r,2) 

8/#  8 2ka>e|s|yo3cos(0-y) 


Equation  (44)  implies  that  the  power  flow  is  always  radial. 

The  direction  of  power  flow  is,  by  definition,  the  direction 
of  the  ray.  Thus  the  rays  corresponding  to  the  surface  waves 
are  straight  lines.  One  can  also  show  that  the  extremals  for 
least  time  in  an  anisotropic,  homogeneous  medium  are  straight 
lines,  so  that  the  rays  not  only  denote  the  direction  of  power 
flow  but  also  are  the  quickest  possible  paths  for  energy  flow. 
Fermat's  principal  is  thereby  verified  for  this  case. 

However,  the  rays  are  not  normal  to  the  wave  fronts. 
The  wave  fronts  are  defined  as  in  the  Isotropic  case  either 
as  (a)  the  contours  of  constant  phase  or  (b)  solutions  to  the 
elkonal  equation.  We  proceed  now  to  find  the  equation  of  the 
wave  fronts  and  to  show  that  these  two  definitions  are  equiva- 
lent in  the  present  case.  For  simplicity  we  suppose  that  the 
phase  of  U(k /*Q»Y)  is  independent  of  y,  as  in  the  case  of  a 
magnetic  dipole  source. 

At  a fixed  Instant  of  time  (or  fixed  value  of/®  ) the 
surface  wave  traveling  in  direction y has  advanced  from  the 
origin  a distance  f I Figure  5 shows,  for  example,  four 

neighboring  phase  fronts.  The  set  of  all  these  phase  fronts 
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fora  the  envelope  of  e curve  C which  is  the  wave  front.  To 
prove  that  C is  a wave  front,  first  repeat  the  process  just 
described  for  a slightly  larger  value  of  />,  say  + . 

There  results  a second  enveloped  curve  . To  show  that  both 
C and  are  wave  fronts,  observe  that  the  distance  between  them 
measured  at  any  point  along  the  normal  to  either  of  them  is 
(correct  to  first  order  in  A/°)  a/7/?0,  where  /d»Q  is  evaluated 
for  the  direction  Y of  this  normal.  This  establishes  the  wave 
fronts  as  a set  of  curves  whose  orthogonal  trajectories  or  wave 
normals  separate  any  two  nearby  curves  by  distances  proportional  to 
the  wave  length  in  the  normal  directions. This  means  that  they 
satisfy  the  eikonal  equation.  We  proceed  to  derive  the  equation 
of  the  wave  fronts: 

The  equation  of  the  phase  front  traveling  in  direction 


(46) 


xcosY  + ysinY  - /*//d 


These  lines  are  the  envelope  of  the  wave  front. 

The  wave  front  is  thus  given  by  the  intersection  of  (46) 


and  its  derivative: 


(47) 


-xsinY  + ycosY 


The  equation  of  the  wave  front  is  therefore 


(48) 


- cos  V + sin  Y 
So  S* 


/6*  cosy  + Siny 
0 


FIGURE  5 

Showing  Wave  Front,  given  In  polar  coordinates  by 
rs*0C09 (0~K)  ■ const.,  as  Envelope  of  Plane  Phase 
Fronts . Each  Phase  Front  Propagates  away  from  the 
Origin  with  its  Characteristic  Phase  Velocity 
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FIGURE  :7 


MIXED  CASE:  X - 1,  XQ  - -4.  Showing  Two 

Wave  Fronts  with  Different  Phase  Velocities 

Y + and  ¥ - . Isotropic  Source  Radiating 
o o 

into  First  Quadrant  Only. 


SQ  + *o  tanV 

The  letter  may  be  rewritten  as 


(51) 


tan(0-Y) 


So 


Comparison  of  (51)  with  (37)  shows  that  0 is  the  same  as  the 
stationary  phase  observation  angle.  Also  (51)  implies  that 

(52)  co»«-y)  - ■/’° 

♦ /»'2 

so  that  (49)  may  be  written  as 

(53)  r^eoad-V)  - J°  . 

Comparison  with  (39)  shows  that  the  wave  front  given  by  (48) 
in  rectangular  coordinates  or  (49)  in  polar  coordinates  is 
also  the  contour  of  constant  phase  in  the  far  field  of  the 
surface  wave.  Definitions  (a)  and  (b)  are  therefore  both 
verified. 


Zt  is  instructive  to  compare  the  propagation  of 
surface  waves  over  an  anisotropic  boundary  with  the  propagation 


I 


of  plane  waves  in  a homogeneous,  anisotropic  medium  such  as  a 
crystalline  medium.  The  latter  is  described  in  considerable 
detail  in  Planck's  "Theory  of  Light"  £3]  . 

One  basic  difference  between  the  two  phenomenon  is 
that  the  propagation  directions  of  the  surface  waves  all  lie 
in  the  plane  of  the  surface,  while  many  of  the  interesting 
features  of  the  crystalline  propagation  (such  as  the  existence 
of  the  primary  and  secondary  optical  axes)  depend  on  the  three 
dimensional  properties  of  the  anisotropic  medium.  If  we  restrict 
the  latter  to  a two  dimensional  medium,  however,  some  interesting 
comparisons  can  be  made.  In  each  case,  a single  infinite  plane 
wave  propagates  energy  in  a direction  which  does  not  coincide 
with  the  normal  to  the  wave  front  unless  the  latter  is  in  a 
principal  direction.  In  each  case,  energy  from  a point  source 
propagates  along  the  radial  lines  from  the  source  although  the 
wave  fronts  are  not  concentric  circles.  However,  the  dependence 
of  phase  velocity  on  direction  is  more  complicated  in  the  surface 
wave  case  than  in  the  crystal  case.  While  in  the 

crystalline  medium  the  wave  fronts  are  ellipses;  on  the  aniso- 
tropic surface  the  wave  fronts  are  more  complicated  curves, 
examples  of  the  latter  are  shown  in  Figures  4a,  b,  c and  7. 

In  the  crystalline  medium,  apart  from  certain  exceptional 
directions,  exactly  two  waves  can  and  do  propagate  from  a 
source.  These  are  cross  polarised  from  each  other,  and  have 
different  phase  velocities.  In  each  wave  E and  H are 
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perpendicular  to  S the  direction  of  energy  flow.  On  the 

other  hand,  on  the  anisotropic  surface,  in  each  direction 

« ; 

either  aero,  one,  or  two  waves  can  propagate  each  with  its 
own  phase  velocity.  Each  wave  is  a hybrid  mode,  a combination 
of  IE  and  TM  modes.  Each  of  the  latter  is  a separate  solution 
to  Maxwell's  equations  but  only  a certain  proper  combination 
of  the  two  satisfies  the  boundary  conditions.  This  same  com* 
bination  guarantees  that  the  surface  wave  energy  flow  from  a 
point  source  will  be  radial  (i.e.,  that  Equation  (43)  holds). 

The  energy  in  the  hybrid  mode  is  not  the  sum  of  the  separate 
energies  in  the  TE  and  TM  parts  so  that  there  is  energy 
coupling  between  the  latter.  E and  H are  not  usually 
perpendicular  to  S. 

The  radiated  field  reflected  from  the  anisotropic 
boundary,  on  the  other  hand,  has  none  of  the  properties  of 
propagation  in  an  anisotropic  medium.  For  the  radiated  field, 

E , H,  and  S are  an  orthogonal  set  and  S is  perpendicular 
to  the  wave  front  which  is  spherical.  While  the  radiated  field 
generally  has  both  TE  and  TM  modes  (even  if  the  source  is  pure 
TE  or  pure  TM),  the  power  flows  of  these  modes  are  orthogonal  to 
each  other. 

EFFICIENCY  IN  A PRINCIPAL  DIRECTION  COMPAQ  «tth  nnjjgN 

When  y ■ 0,  then  from  (17),  y^  ■ 0 and,  from  (37), 

0 » 0.  We  observe  from  (16)  that  YQ  • X or  -1/XQ  (if  these 
are  positive).  Substitution  of  -1/XQ  for  YQ  and  zero  for  V 
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into  (35)  yields  ^ ■ 0.  Hence  there  is  no  excitation  o£ 
the  mode  for  which  » -1/X0  the  direction  Y « 0 even  if 
this  mode  appears  to  be  a free  solution  of  Maxwell's  equations 
and  the  boundary  condition.  For  the  YQ  « X mode  at  V'  “0, 


(54) 


2kX* 

v/l+X2 


From  (17),  correct  to  second  order  in  Y 


(55) 


Vo  + 


V2(xo-x)(hx2) 

1+XX_ 


So  that  from  (37),  correct  to  first  order  in  Y 


(56) 


From  (38) 


2X(X-X  ) 

0 - Y (1  + — 

1+XX_ 


(57) 


a(0)  - - Tu?  (1  + 


2X(X-X  ) 


Finally  from  (45) 

_ 2,kX|U(k/fo>0)|2 


(58) 


C0€  (1+ 


For  the  purpose  of  comparing  our  result  with  Cullen's  [l] 
let  us  define  the  efficiency  4 of  surface  wave  excitation  in  a 
given  direction  0 as  the  ratio  of  surface  wave  power  radiated 
outward  between  the  meridian  0 and  0+d0  to  the  total  power 
radiated  outward  between  these  meridians  in  the  limit  as  d0 
approaches  zero.  From  (31)  and  (57),  we  see  that  for  0 ■ 0 


j-j-w  Vi*** 


i&iXafca&i  - a tsI  i^-  - ^r^^^e^paigiisaisfea^a^wisaftiaa^^  ~.ia|£4»i^:  7;:.  :;*$&>»«■»  ij| 


(59) 


where 

(60) 


P./P. 


+ Ps/Pr 


*8  4yx/u(k^Q>0)/2  

^r  2X(X-X  ) i 2 ^ 

(1+ 2—)  ( sin0d0  / F(0,O)+R1  F(ir-0,O)f 

1+xx°  1 

where 

2X(iX  cosQ-l) 

(61)  R,  - 1 + 8 • 

A XQcos*0  - icos0(XXQ-l)+X 

Apart  from  the  fact  that  Cullen's  formula  for  efficiency 
applies  only  to  the  case  XQ  » 0 (single  corrugations),  there 
are  three  differences  between  our  formula  and  his. 

(1)  A factor  c.  1+X  is  missing  from  the  denominator  of 
(60) . The  absence  of  this  factor  appears  to  be  an  inherent 
difference  between  the  two  and  three  dimensional  problems,  since 
its  counterpart  is  also  missing  in  the  rigorous  expression 
(130)  [ 2 } for  the  surface  wave  in  the  three  dimensional  iso- 
tropic case.  The  difference  between  the  two  and  three  dimen- 
sional problems  in  this  respect  may  be  traced  back  to  the 
difference  between  the  three  dimensional  integral  representa- 
tion for  the  field  as  given  by  (22)  and  Cullen's  integral 
representation  (19)  of  [lj.  The  expressions  are  almost  the 
same,  but  the  factor  p’  in  the  integrand  of  (22)  has  no 
counterpart  in  (19)  of  flj.  This  factor  p is  the  Jacobean 
of  the  change  of  coordinates  in  the  transform  plane  from 


-v  ^ jjJ-,  .is- , 


rectangular,  with  area  element  dudv,  to  polar,  with  area  element 
pdpd y . (This  change  of  coordinates  is  carried  out  explicitly 
in  £ 2 ] ).  The  same  change  of  coordinates  in  the  two  dimensional 
case  need  never  be  explicitly  considered  since  it  degenerates 
to  the  identity  transformation. 

When  the  residue  theorem  is  applied,  p is  replaced 
by  which^at  0 « 0? equals  Jl  +X^  . When  the  latter  factor 

is  squared  for  the  computation  of  power,  the  factor  of  lfX  , 
which  appears  in  the  numerator  of  Cullen's  expression, (and 
would  otherwise  also  appear  in  (60)) is  cancelled  (Note  that 
Cullen  uses  j*  for  our  X) . 

(2)  In  computing  total  radiated  power,  we  must  integrate 
over  a large  sphere,  not  over  a cylinder  as  in  Cullen's  case, 
so  that  the  factor  sin0  appears  in  the  integrand  in  (60)  but 
does  not  appear  in  the  analogous  expression  (35)  of  [l]  . 

Again  this-  factor  is  a Jacobean  which  is  not  the  same  in  two 


and  three  dimensions. 
(3)  The  factor  1 + 


2X(X-X  ) 


appearing  in  the  denominator 


of  (60)  does  not  appear  in  Cullen's  expression.  It  arises 
from  the  fact  that  the  wave  number  k/TQ  is  a function  of 
direction  V , i.e.,  that  the  boundary  is  anisotropic,  so 
that  s is  a function  of  0.  In  Cullen's  case  s equals 
m /9q  and  is  independent  of  0.  To  say  the  same  thing  in  more 
physical  though  perhaps  imprecise  terms:  the  fact  that  the 
wave  fronts  of  a point  source  are  not  circles  implies  that. 


35. 


in  the  stationery  phase  evaluation,  the  far  field  will  vary 
with  the  curvature  of  the  wave  front  at  the  stationary  phase 
point  - the  flatter  the  phase  front  and  the  mere  "stationary" 
the  stationary  phase  point,  the  larger  the  contribution  to  the 
far  field.  This  phenomenon  has  no  analogue  in  the  isotropic 
boundary  case. 


Proof  of  (43): 


Ie  Is  sufficient  to  prove 


(Al) 


ten(|- Y) 


21  yprir2 


:1  “ r2 


or  in  view  of  (37) 


(A2) 


i+ 


21rlr2 


rl  “r2 


From  (35) 

(A3)  2ir,r0  2SC(X-XJ  y (XC2-X  S2+  / XXJ 


rl  “r2 


(xc  +XQS  + yoxxo)‘  + s‘c‘yo"(x-xo) 


where  we  abbreviate 


(AA) 


S-sinV, 


C » cos  Y 


Since  Yq  is  the  root  of  the  denominator  in  (11)  we 


differentiate  this  denominator  to  show  that 
(A5)  ^(2  yQ(X0C2+XS2)  + l-XXft)+2(l+  /2)(X-Xn)SC  - 0 . 


This  implies  that 


(A6)  - V • 


2(X-XJSC 


2(X-X  >SC 


C 2yo<XoC  + “*>  + t/(l+XX0)2+(X-X0)2sin22Y 

where  the  sign  alternative  is  to  be  chosen  to  insure  equality. 

We  find  then  that  (A2)  will  hold  if  and  only  if 

(a?)  yo(xc2+x0s2+  jroxxo)(2>'o<xoc2+xs2)  + i-xx0)  - 
(xe2+x„s2+  y„xx„)2  + s2c2  y„2(x-x„)2  , 


4 


) 


! 


i 

i 

\ 
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(A7)  may  be  written  as 
(A8) 

y03(2XX0)(X0C2+XS2)+  yo2(p0-2(xxo)2  + 2(xoc2+xs2)(xc2+xos2) 
-s2c2(x  -x)2)  + y0(xc2+x0s2)(i-3xx0)  - (xe2+x0s2)2  . o . 

Since  from  (16) 

(A9)  yo2(xQc2+xs2)  - xc2+x0s2  + *0(xxo-i)  9 

we  may  substitute  for  the  cubic  term  in  )fQ  in  (A8)  reducing  this 
to  terms  of  second  degree  or  less,  which  upon  combination  with 
the  remaining  terms  of  (A8)  yield 

(A10)  yo2((X02+X2)S2C2+XX0(2C4+2S4+2S2C2-l)  ) + (1-XXQ)  yo(XC2+XoS2) 
-(XC2  + XQS2)2  - 0 

The  coefficient  of  )^2  may  be  factored  by  (XC2+XQS2),  and  with 
the  cancellation  of  the  latter,  (A10)  reduces  to  (A9) . Since 
all  the  algebra  is  reversible,  we  have  proved  (43). 


PART  II 


39. 


EXTENSION  OF  SOME  RESULTS  OF  CULLEN  ON 
THE  EXCITATION  OF  PLANE  SURFACE  WAVES" 


In  [l}  , Cullen  has  obtained  a formula  for  the  excita- 
tion efficiency  of  surface  waves  launched  over  a metal  clad 
dielectric.  The  formula  is  based  on  the  following  assumptions. 

(a)  the  dielectric  thickness  d is  small  compared  to 
the  wavelength  X 

(b)  the  energy  carried  by  the  surface  wave  in  the  dielectric 
is  negligible  compared  to  the  energy  carried  by  the  surface  wave 
in  the  space  above, 

(c)  an  infinite  horizontal  magnetic  line  current  is  the 
source. 

For  one  application  encountered  at  Technical  Research 
Group  it  was  of  interest  to  obtain  a formula  generalizing 
Cullen' s which  drops  the  first  restriction  and  replaces  his 
source  by  a finite  horizontal  magnetic  line  current.  A restriction 
on  thickness  is  still  invoked  namely 

(1)  d4X/(2V?T) 

but  this  is  used  to  insure  that  only  a single  surface  wave  can 
propagate.  When  the  thickness  is  not  small  compared  to  the  wave- 
length, the  assumption  that  the  energy  carried  by  the  surface 
wave  in  the  dielectric  is  negligible  is  no  longer  tenable  and  it 
is  dropped  here.  It  is  rather  surprising  to  discover  that  despite 
the  apparent  immediate  increase  of  complexity  caused  by  the  waiver 
of  these  assumptions,  our  final  formula  for  excitation  efficiency 


it  hardly  more  complicated  than  Cullen's. 


we  simply  point  out,  line  by  line,  what  changes  must  be  made  in 
Cullen '8  derivation  when  assumptions  (a)  and  (b)  are  not  made,  and 
as  a first  step,  we  obtain  an  Intermediate  formula  when  only  these 
assumptions  are  dropped.  In  order  to  change  assumption  (c),  we 
then  rely  on  the  results  of  the  three  dimensional  analysis  given 
in  reference  [2]  and  in  the  first  part  of  this  paper. 

We  use  on  the  whole  Cullen's  notation.  The  first  ex- 
ception is  that  we  let  € be  the  relative  dielectric  constant 
of  the  dielectric.  However,  in  contrast  to  the  first  part  of 
this  paper,  now  the  implicit  time  dependence  is  €rwt  and  0 is 
the  grazing  angle  rather  than  the  incidence  angle. 

REMOVAL  OF  ASSUMPTIONS  (a)  and  (b) 

The  reflection  coefficient  for  a TM  plane  wave  grazing 

the  metal  clad  dielectric  surface  at  angle  0 is 
g . isinO  + f(cos9) 
jslnO  - X(cosO) 

where 

(2)  f(t)  - tan  (kd/e*tz) 

Section  (10)  of  Cullen  may  now  be  read  through  expression 
(44)  if  the  latter  is  replaced  by  kX(cos9).  Cullen's  formula 
for  the  radiated  power  (34)  remains  correct  with  his  f replaced 
by  X(cosO). 

In  order  to  find  the  formula  for  the  surface  wave  power, 
observe  that  equation  (19)  of  £l^  remains  correct  if  k is 


replaced  by  kf(^/k).  The  generalized  reflection  coefficient 

«>  «j»>  - SSfgfB- 


then  appears  in  (19).  Its  residue  at  the  pole  u ■ kp  where  /Q* 

A -J>W  ■ 1” 

(4>  frk^oA-zfr2) 

^(€-1)  [ek^+k^k^s+D+k^J 

The  leading  term  of  (20)  [l]  is  then 

(5)  ^ p 

(«-l)[  €k2  + kjdtk^Cc+D+k2)  7 

The  power  radiated  in  the  surface  wave,  is  then  (instead  of  (36) [l]) 

(6)  */?  A'iklh.\V«v1¥ 


k «+kld(k1  (e+l)+k  ) 


dy  - l/2k,  . 


cosyy  k €-/ 


2ccos(d  / k e-  / d L *) 


1 
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Let  us  now  Introduce  a simplifying  notation 
(9)  * icn»  “ k/»*-l  ) 


ki  - ky/n  -1 

where  n is  the  equivalent  "index  of  refraction”.  The  surface 
wave  existence  condition  (43)  of  £lj  may  be  written  as 

(10)  e ynZ-l  » >/«-n^tan(kdy€-n^)  , 
so  that  (8)  becomes 

(11)  i . d(c»l)  (n^(c-fl)-c)  + U/k)Jnz-l  . 

2 2e(e-n2) 

The  ratio  of  the  power  in  the  surface  wave  carried  above  the 
dielectric  to  the  power  carried  in  the  dielectric  is 


(12) 


iig-n2). 


2 ,/n2-l  (£d(€-l)(n2(€+l)-€)  + ej n^-1  ) 

When  kd<£^l  so  that  n ^ 1 this  ratio  is  large,  and  the 
term  I2  may  be  neglected  in  (6)  compared  to  1^.  In  the  general 
case,  however,  I2  is  not  negligible,  and  (6)  must  be  written 

as 

2k V2e~2kh^Z_1  l(e-n2) 

n(e-l)  ^€+(kd/nz-l)  (n2(c+l)-e)J 
From  (35)  of  [l]  the  radiated  power  is 
(H) 


Pr  - 1^7  *0  V2 


»/2 

5 


F2(0)de 


and  the  efficiency  is 


(15) 


? ■_,r+Ts 


where 

(16) 


-2kh 


,i 


n<€-l)^€+(kdJ?^l)(n2(e+l)-€)]^2F2(e)d0 


FINII 


,INE  SOURCE 


Two  changes  must  be  made  in  formula  (16)  if  the  line 

source  is  finite  instead  of  infinite.  Consider  any  source 

whose  three  dimensional  far  field  pattern  in  free  space  is 

uniform  In  the  xy  plane,  i.e.,  independent  of  angle,  and  has 

H perpendicular  to  this  plane.  The  analysis  of  [2]  - or  of 

the  proceeding  section  - shows  that  an  extra  factor  of  n appears 

2 

in  the  surface  wave  field  and  hence  n in  the  surface  wave 
power.  In  the  radiation  field  the  power  flow  must  be  integrated 
over  a large  sphere  whose  area  element  is  ordinarily  sin9dOd0  - 
but  since  9 is  now  the  grazing  angle  instead  of  the  Incidence 
angle  in  the  present  notation  - this  factor  is  cos9d9d0  here. 

The  d0  cancels  the  same  factor  in  the  surface  wave  power  ex* 
pression,  if  we  consider  the  ratio  of  power  flows  per  unit 
azimuthal  angle.  The  correct  expression  for  efficiency  is  thus 
(15)  with  (16)  replaced  by 
(17)  m — 


Is  m.IXl  f.-A.-a*#1* 

(€-1)  £ €+(kd/n2-l)(n2(€+l)-€)J  j'1 


cos9F*(9)d9 


Both  additional  factors  in  the  three  dimensional  case  tend  to 


enhance  the  power  in  the  surface  wave  at  the  expense  of  power 
in  the  radiation  field  when  compared  to  the  two  dimensional 
case. 

It  is  also  well  to  point  out  that  the  preceeding  formulas 
for  efficiency  will  in  general  be  gravely  in  error  if  the  source 
is  not  symmetric  with  respect  to  a plane  parallel  to  the  ground 
plane,  and  can  be  considerably  in  error  if  the  source  is  anything 
but  isotropic  in  the  xy  plane.  This  is  especially  true  for  those 
high  values  of  efficiency  which  are  primarily  due  to  the  cancela- 
tion of  the  directed  and  reflected  rays  in  the  radiation  field. 
(See  Fig.  10  of  fl]  ). 

Figure  9 shows  the  two  and  three  dimensional 

efficiencies  for  four  values  of  e vs  h/A  and  d/x,  the  height 

of  the  source  and  the  thickness  of  the  dielectric,  respectively, 

2 

in  wavelengths.  The  corresponding  values  of  n are  shown.  The 

values  of  h/X  where  the  radiated  field  is  zero  at  both  45° 

and  60°  grazing  angle  are  shown.  These  are  not  necessarily  the 

most  efficient  heights.  If  the  conclusions  of  reference  [5}  were 

to  be  taken  literally  for  the  dielectric  clad  ground  plane  then 

at  the  former  height  the  efficiency  would  always  be  100%.  Cullen's 

2 

conclusions  based  on  the  case  e - 2.56,  n >1.25  also  appear 
not  to  be  highly  relevant  for  other  values.  The  optimum  height 
for  maximum  efficiency  is  apparently  often  zero.  Moreover, 
the  optimum  height,  which  is  always  a small  fraction  of  a wave- 
length, is  not  critically  different  from  zero.  The  efficiency 
does,  however,  depend  critically  on  d/x  (or  n).  None  of  these 
conclusions  are  apparent  from  Cullen's  work. 


:ure  9c 


efficiency  * of  Excitation  of  Surface  Wave  over 
Dielectric  Slab  of  Thickness  d/x  for  Magnetic 
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ERRATA  SHEET 

Page  33*  Lina  9 and  Page  34*  Line  8*  replace 
"1  + X2”  by  " Jl  + X2  " 

Page  33*  Line  7*  ac  end  of  sentence  insert: 

"The  stationary  phase  evaluation  of  the  V Integration 
reduces  this  factor  to  (1+X2)*^.  " 

Page  43,  Line  9 - replace  "n"  by  "\/n"”  " 

Page  43,  Line  10  - replace  "n2"  by  "n" 

Equation  (17)  - delete  factor  "n"  in  numerator. 

To  the  references  add: 

[7]  Fernando  and  Barlow*  "An  Investigation  of  the 
Properties  of  Radial  Cylindrical  Surface  Waves 
Launched  over  a Flat  Reactive  Surface", 

Proc.  of  the  Inst,  of  Elect.  Eng.*  Part  B* 

May  1956*  p.  307. 
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